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The out-of-equilibrium transport properties of a double 
quantum dot system in the Kondo regime are studied theo- 
retically by means of a two-impurity Anderson hamiltonian 
with interimpurity hopping. The hamiltonian is solved by 
means of a non-equilibrium generalization of the slave-boson 
mean-field theory. It is demonstrated that measurements of 
the differential conductance dl /dV , for the appropriate values 
of voltages and tunneling couplings, can give a direct obser- 
vation of the coherent superposition between the many-body 
Kondo states of each dot. For large voltages and arbitrar- 
ily large interdot tunneling, there is a critical voltage above 
which the physical behaviour of the system again resembles 
that of two decoupled quantum dots. 



Recent experiments |l] |3|] have shown that new physics 
emerge when the transport properties of quantum dots 
(QD's) at temperatures (T) below the Kondo tempera- 
ture (Tk) are studied. Q QD's offer the intriguing pos- 
sibility of a continuous tuning of the relevant parameters 
governing the Kondo effect [|| as well as the possibility 
of studying Kondo physics when the system is driven out 
of equilibrium in different ways. || These experimental 
breakthroughs have opened up a new way for the study 
of strongly correlated electrons in artificial systems. The 
Kondo anomaly appearing in the density of states (DOS) 
of the QD reflects the formation of a quantum-coherent 
many-body state. Motivated by the recent experimental 
advances in the study of double quantum dots (DQD) 
M it is thus interesting to study what happens when two 
QD's in the Kondo regime are coupled. Previous theoret- 
ical studies of this problem at equilibrium have focused 
on the competition between Kondo effect and antiferro- 
magnetic coupling generated via superexchange [pL[| or 
via capacitive coupling between dots. 
In this Letter we focus on the study of a DQD in the 
Kondo regime driven out of equilibrium by means of a 
DC voltage bias. There have hitherto been only few at- 
temps to study this problem |11| but a clear picture of 
the problem is yet missing. Following the recent work of 
Aono et al [[l2| and Georges and Meir [|| we employ the 
slave boson (SB) technique [ p"3| in a mean field approx- 
imation (MFA) and generalize it to a non-equilibrium 
situation. This MFA allows us to include nonperturba- 
tively the interdot tunneling term (i.e coherence between 
dots). The different physical regimes that appear as the 
ratio r c = £c/r changes (tc is the interdot tunneling 
coupling and T is the single particle broadening com- 
ing from the coupling to leads (IJ]) can be explored by 



measuring the non-linear transport properties of the sys- 
tem. Our results can be summarized in Figs. 1 and 2: 
the differential conductance dl/dV of the DQD directly 
measures the transition (as r c increases) from two iso- 
lated Kondo impurities to a coherent superposition of the 
many-body Kondo states of each dot, which form bond- 
ing and anti-bonding combinations. This coherent state 
which occurs for r c > I is reflected as a splitting of the 
zero-bias anomaly in the dl/dV curves. This splitting de- 
pends non-trivially on the voltage and on the many-body 
parameters of the problem. For large voltages, we find 
that there is a critical voltage above which the coherent 
configuration is unstable and the physical behaviour of 
the system again resembles that of two decoupled QD's, 
i.e two Kondo singularities pinned at each chemical po- 
tential, even for r c > 1. This instability is reflected as a 
drastic drop of the current leading to singular regions of 
negative differential conductance (NDC). 
Model: In typical experiments, £/i ntr adot, Ae >> T 
(£Antradot is the strong on-site Coulomb interaction on 
each dot, Ae is the average level separation), which al- 
lows one to consider a single state in each QD. [[l5| We 
can model the DQD with a (N=2) fold degenerate two- 
impurity Anderson hamiltonian with an extra term ac- 
counting for interdot tunneling. Each impurity is coupled 
to a diferent Fermi sea of chemical potential fj,r, and [1r, 
respectively. In the limit [A n t ra dot — > oo (on each QD) 
and C/i n tcrdot - ► flql , the hamiltonian may be written 
in terms of auxiliary SB operators |r| plus constraints: 

+ ^ E(/i + fhMb{f La ) 

-J= Y V <*( C k a ,a b tf<*a + fL b a C ka ,a) 

X a (Y^J a , + bib a -l). (1) 
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c\ a (ck a .cr) are the creation (annihilation) operators for 
electrons in the lead a. To simplify the notation we con- 
sider henceforth that Vl = Vr = Vq and e£ CT = e^ CT = £q 
(i.e, Tk is the same for both dots at equilibrium. The 
generalization to different Tk 's is straightforward). The 
even-odd symmetry is broken by the interdot coupling 
tc- In the SB representation, the annihilation opera- 
tor for electrons in the QD's, c aa is decomposed into 
the SB operator b' a which creates an empty state and a 
pseudo fermion operator f a<7 which annihilates the singly 



occupied state with spin a in the dot a: c aa — > b^ a f aa 
( c aa ~~ > faaba)- In the last term of (1), the charge oper- 
ator Q a = J2 a faafaa + b\b a has been introduced. This 
term represents the constraint Q a = 1 in each dot with 
Lagrange multiplier A Q . This constraint prevents double 
occupancy in the limit U — ► oo. 

Solution: In the lowest order, we assume that the 
SB operator is a constant c- number b a {t)/\N — 
(b a (t))/VN = b a neglecting the fluctuations around the 
average (b a (t)) of the SB. At T=0, this MFA is correct 
for describing spin fluctuations (Kondo regime). Mixed- 
Valence behavior (characterized by strong charge fluc- 
tuations) cannot be described by the MFA. This re- 
stricts our non-equilibrium calculation to low voltages 
V « e . Charge fluctuations can be included as ther- 
mal or quantum fluctuations (l/N corrections). [ ^3|Jl7| 
Defining V a = Vob a and tc — ic&L&R we obtain from 
the constraints and the equation of motion of the SB op- 
erators the selfconsistent set of four equations with four 
unknowns (b~ L ,b R , X L , X R ): 

b 2 L(R) + ^2(fl(R)<rfL(R)<r) = TV 

(T 

Yi<B r /J f \ 

+ ~§ ^l^R(L)afL(R)a) + ^L(R)b\ {R) = (2) 
cr 

In order to solve (2) we need to calculate the non- 
equilibrium distribution functions: G^ak ,a(t~t') = 

i(4 a , CT (*')/a<r(t)> and G< aa ,Jt-t>) EE i(fl, a W)faa(t)). 

They can be derived by applying the analytic contin- 
uation rules of Ref. |l^] to the equation of motion of 
the time-ordered Green's function along a complex con- 
tour (Keldysh, Kadanoff-Baym or a more general choice) . 
This allows us to relate G< h Jt-t') with G < , (t - 

t<) and Gr ata , a (t - t>) ee -iO(t - t')({f aa (t)jl >a (t')}) 
and close the set of equations (2) in Fourier space: 

^HR) _ ■ f fe r < , x _ 1 

T J 2ir e > ~ N 

-^P-(h{R)-e )=iJ ^G< L{RR) {e){e-e L{R) ), (3) 

with e a = eo + A Q and r a = b^T. For tc = , e Q and 
r a give, respectively, the position and the width of the 
Kondo peaks in the dot a (at equilibrium, and in the 

Kondo regime, ^Ji 2 a + T 2 a = T£ = De^ e °V r ). g The 
distribution functions in the QD's are: G^ L , RR ^(e) — 

2*(f 'l(r) fun) (z)[(£-ZR(L)) 2 +r 2 R(L) ]+f R{ L)fR(L)(t)i 2 c ) 
[(e-lL+iT L )(<L-t R +iT R )-t 2 c ][(e-e L -iT L ){t-t R -iT R )-t 2 c ] ' 

/l(_r)( £ ) is the Fermi-Dirac function in the left (right) 



lead. Note that the presence of t c in the denomina- 
tors indicates that the interdot tunneling enters nonper- 
turbatively in the calculations and, then, coherent ef- 
fects between dots are fully included. Due to the in- 
terdot tunneling, the Kondo singularities of each dot 
at £l and e R combine into coherent superpositions at 

e± = \{{tL + cr) ± \J(cl - £i?) 2 +4^}. Of course, at 
equilibrium b^ — b R — b, Xl = X R = A, we recover the 
results of Refs. p2|,p[' Note that the formation of coher- 
ent superpositions of the Kondo singularity is not triv- 
ially related with its single-particle counterpart (forma- 
tion of bonding and antibonding states at eo ± tc)- Let's 
focus for simplicity in the equilibrium case (ej, = e R ), 
the splitting is given by 5 = e + — e_ = 2tc which is 
a many-body parameter (given by the strong renormal- 
ization of the interdot tunneling due to the Kondo ef- 
fect). S depends non-linearly on the single-particle split- 
ting So — 2tc (see Inset of Fig. 3a). In the Kondo limit, 
{[(e + t c ) 2 + f 2 ][(e - t c ) 2 + f 2 ]}V4 = T^e^r-i). 
From the solution of Eq. 3 we obtain the current 
I = !ME* tl , V L G<^ kLa {t, t)} and DOS in each QD: 

Results: We solve numerically (for T=0) the set of non- 
linear equations (3) for different voltages = V/2 and 
Hr = —V/2, eo = —3.5,1? = 60 (Kondo regime with 
Tj( ~ 10 -3 ) and different values for the rest of parame- 
ters (all energies in units of T) . Depending on the ratio 
t c = tc /r we find two different physical scenarios for 
r c < 1 and t c > 1. 
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FIG. 1. I-V curves for different values of r c < 1 and 
eo = —3.5. Caption: dl/dV curves for the same parameters. 

In Fig. 1 we plot the I-V curves (for clarity, we show 
only the V > region) for r c < 1. The two main fea- 
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tures of these curves are: i) An increase of the linear 
conductance Q = dl /dV\v=o as t c increases; ii) a satura- 
tion, followed by a drop, of the current for large voltages. 
This drop sharpens as t c -> 1. These features are more 
pronounced in a plot of the dl/dV (inset of Fig. 1). 
As t c increases, the zero-bias anomaly (originating from 
the Kondo resonance in the DOS of the dots) becomes 
broader and broader until it saturates into a flat region 
of value 2e 2 /h (unitary limit) for r c = 1. The reduction 
of the current at larger V is reflected as NDC regions in 
the dl/dV curves. For r c = 1 this NDC becomes singular. 
For t c > 1, and contrary to the previous case, Q decreases 
for increasing values of r c (Fig. 2a). This reduction of Q 
can be atributed to the formation of the coherent super- 
position of the Kondo states. This can be clearly seen 
as a splitting A = 25 in the dl/dV curves (Fig. 2c): In- 
creasing t c , the zero-bias conductance decreases whereas 
two maxima at ±V pea k show up (the arrow shows the 
splitting A = 2V pea k for the maximum value of r c in the 
figure). Fig. 2c demonstrates that the dl/dV curves of 
a DQD in the Kondo regime directly measure the coher- 
ent combination between the two many-body states in 
the QD's. For larger voltages, the sharp drop of the cur- 
rent (Fig. 2a) reflects as strong NDC singularities in the 
dl/dV curves (Fig. 2b). 
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FIG. 2. a) I-V curves for different values of r c > 1 and 
eo = —3.5. b) dl/dV curves for the same parameters, c) 
Blow up of Fig. 2b. The arrow shows the splitting A = 28 
for Tc — 1.5. 

The position of these singularities moves towards 
higher |V| as t c increases. In order to explain the re- 
sults of Figs. 1 and 2, we plot in Fig. 3a e± as a function 
of V > for different values of t c . For r c = (thick 
solid line) , this corresponds to a plot of and e R (i.e 
the positions of the Kondo resonances for the decoupled 
QD's) as a function of V. We obtain, as expected, that 
each Kondo resonance is pinned at the chemical potential 
of its own lead, cl = = V/2 and e R = [1r = —V/2. 



As the interdot coupling is turned on, the voltage depen- 
dence becomes strongly non-linear. For low V, the curves 
for r c 7^ do not coincide with the curves for t c = (i.e, 
I^l/r)- This situation, however, changes as we increase 
V; the level positions e± converge towards the chemical 
potentials \i L / R in a non-trivial way. 
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FIG. 3. a) e± vs. V for different values of r c (eo = —3.5). 
Inset: Many body splitting (8) as a function of the single par- 
ticle splitting (<5 ) for V = and e = -2.0, -2.5, -3.0, -3.5 
(from top to bottom), b) Fl vs. V. c) Fr vs. V. d) tc vs. V. 

The voltage V for which e+ — e~ coincides with the 
chemical potential difference V gives the position of the 
peak in the positive side of the dl/dV (Fig. 2c). This 
voltage is the solution of the equation 5(V pea k) — V^ oa k 

where 5(V) = yj (e L - e R ) 2 + U 2 C , with e L / R given by 
Eq. (3). Note the implicit (and non-trivial) voltage de- 
pendence of S(V). Tl, Fr and tc follow a similar be- 
havior as a function of V (Figs. 3b, 3c and 3d). For 
V > Vp ea k, we find numerically that S(V) w V, a rela- 
tionship that becomes asymptotically exact as V — > oo. 
The equation S(V) — V, has stable solutions tc f° r 
[ (ei,-g R ) j2 < 1; whilc for [ (^- g «) ]2 > 1, the only stable 
solution is tc = 0, corresponding to current 1 = 0. Wc 
denote the crossover voltage where [ ( ^~ £ ' r) ] 2 = 1 by V*. 
For finite voltages V > Fp 0a k, on the other hand, the 
relation 6(V) = V is only approximate, so that at the 
crossover ss V* , the quantity tc and hence / drop to a 
much smaller, but still finite, value instead. Nevertheless 
the crossover at V ~ V* still indicates the beginning of 
the NDC region. 

To illustrate this, we plot in Fig. 4 the left and right 
QD's DOS for r c = 1. At equilibrium (V=0), the Kondo 
singularity at e = splits into the e± combinations. For 
V/Tft = 2 the coherence is still preserved but the physi- 
cal picture utterly changes for higher voltages (V/T^ = 4 
and V/1% = 6). In this case, the previous configuration 



3 



is no longer stable, the coherence between dots is lost 
(tc — > 0), the dots are almost decoupled and the Kondo 
resonances in each dot are pinned again at their own 
chemical potential: the weight of the left (right) DOS at 
e = Hb.il) is almost zero (even though r c = 1). 
This instability resembles that of the SB at T ^ in 
the single impurity Anderson hamiltonian. 1 17| 1^] In the 
MFA the SB behaves as the order parameter associated 
with the conservation of Q. When b ^ the gauge sym- 
metry b — > be 10 , / — > fe l6) associated with charge conser- 
vation is broken and the MFA has two phases 6^0 and 
b = separated by a second order phase transition. It is 
important to point out that the fluctuations do not de- 
stroy completely this b — ► behavior (the SB fluctuations 
develop power law behavior replacing the transition by a 
smooth crossover). We speculate that in our problem this 
zero-temperature transition at finite V may be also ro- 
bust against fluctuations but I /N corrections are needed 
to substantiate this argument. Work in this direction is 
in progress. 




FIG. 4. DOS for the left and right dot for r c = 1 and differ- 
ent voltages (eo = —3.5). Caption: Transmission probability 
of the DQD for the same parameters. 

In closing, we have demonstrated that the non-linear 
transport properties (dl/dV) of a DQD in the Kondo 
regime directly measures the transition (as tc increases) 
from two isolated Kondo impurities to a coherent bond- 
ing and antibonding superposition of the many-body 
Kondo states of each dot. While for tc < T the con- 
ductance maximum is at V = 0, for tc > T the trans- 
port is optimized for a finite V matching the splitting 
between these two bonding and antibonding states. For 
large voltages (and tc > T) there is a critical voltage 
above which the coherent superposition is unstable and 
the physical behavior of the system again resembles that 
of two decoupled QD's. This leads a strong reduction of 
the current and singular regions of negative differential 
conductance. Concerning the observability of these ef- 
fects: In our MFA the maximum value of S ranges from 
S ~ 20T£ - 500T£ (inset of Fig. 3a) giving, for the ex- 



periment of Ref. D (r ~ ISOfieV), S ~ 3/xeV - 75/zeV 
(30mK-750mK) which is within the resolution limits of 
present day techniques. 
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